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Abstract 

Processes of heavy quark photoproduction at HERA energies and be- 
yond are investigated using the semihard {k± factorization) approach. 
Virtuahty and longitudinal polarization of gluons in the photon - gluon 
subprocess as well as the saturation effects in the gluon distribution func- 
tion at small x have been taken into account. The total cross sections, 
rapidity and p± distributions of the charm and beauty quark photopro- 
duction have been calculated. The results are compared with ZEUS ex- 
perimental data for charm photoproduction cross section. 
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1 Introduction 



The heavy quark electro- and photoproduction has become an increasingly important sub- 
ject of study 0, H, H]. In particular it originates from come into operation of HERA ep- 
collider energy of which is high enough to achieve the small x physics region. Indeed first 
experimental results obtained at HERA Q| showed new physical phenomena in 7p- and 
ep-interactions at HERA energies; resolved photon interactions, jets in photoproduction 
produced in 'hard parton collisions, significant fraction of events with a large rapidity gap. 
The measurement of the total charm photoproduction cross section by ZEUS collaboration 
Q is a very important stage of researches. 

The heavy quark production at HERA is a very interesting and important subject of 
study. As it is dominated by photon- giuon fusion subprocess (Fig. 1) one can study the 
gluon distribution functions G{x,Q^) in small x region (roug hly at X > 10"''). Secondly 
this last issue is important for physics at future colliders (such as LHC): many processes at 
these colliders will be determined by small x gluon distributions. 

At HERA energies and beyond heavy quark production processes are of the so-called 
semihard type In such processes by definition a hard scattering scale Q (or heavy 

quark mass M) is large as compared to the Aqcd parameter, but Q is much less than the 
total center of mass energies: Aqcd ^ Q ^ \/s. Last condition implies that the processes 
occur in small x region: x ~ M^/s ^1. In such a case the perturbative QCD expansion 
has large coefficients 0(ln"'(s/M^)) ~ 0(ln"(l/x)) besides the usual renormalization group 
ones which are 0(ln"((5^/A^)). So, in perturbative QCD the heavy quark photoproduction 
cross section (as result of photon-gluon fusion process) has the form Q : a^g — a°g + oisf^l^g 
-t- where al^g was calculated by Ellis and Nason The photon - gluon fusion cross 
section in low order decreases vs. s at s — > oo : cr°g ^ A-P / sh\{s / Q'^) . But in the same 
limit al^g — > Const, because heavy quark photoproduction cross section is dominated by 
the contribution of the gluon exchange in the t - channel. These results in breakdown of 
standard perturbative QCD expansion and the problem of summing up all contributions of 
the order (a,, In(QVA^))". ln(l/a;))" and (a^ ln(gVA2) ln(l/a:)" in perturbative QCD 
appears in calculation of a~fg. 

It is known that summing up the terms of the order {as ln(Q^/Af^))" in leading loga- 
rithm approximation (LLA) of perturbative QCD leads to the linear DGLAP evolution equa- 
tion for deep inelastic structure function . Resummation of the large contributions of the 
order of {as ln(l/a;))" leads to the BFKL evolution equation [|l^ and its solution gives the 
drastical increase of the gluon distribution: xG{x,Q^) x^'^°,ujq = (4A'c In 2/7r)as((5§). 
The sharp growth of the parton density as x — > makes the parton-parton interactions very 
important, which in turn makes the so-called GLR evolution equation essentially nonlinear 

The growth of the parton density at a; — > and interactions between partons induce 
substantial screening (shadowing) corrections which restore the unitarity constrains for 
deep inelastic structure functions (in particular for a gluon distribution) in small x region 
pT[ |. These facts break the assumption of the standard parton model (SPM) about x and 
transverse momentum factorization for a parton distribution functions. We should deal 
with the transverse momentum factorization {kj_ factorization) theory 13 1 or with the 



2 



so-called semihard approach |ll| in case of accounting for shadowing corrections at small 
X besides taking into account the virtuality and the longitudinal polarizations of initial 
gluons. 

In semihard approach |ll] screening corrections stop the growth of the gluon distri- 
bution at a; ^ 0. This effect was interpreted as saturation of the parton (gluon) density: 
gluon distribution function xG{x,Q'^) becomes proportional to Q^R^ at < qQ{x) and 
cr ~ ■^xG{x, Q^) ~ R^. The parameters R and qQ{x) can be considered as new phenomeno- 
logical parameters: the R is related to the size of hadron or of black spots in a hadron and 
the parameter is a typical transverse momentum of partons in the parton cascade of 
a hadron, which leads to natural infrared cutoff in semihard processes. The parameter Qq 
increases with s ||ll|, Q . So because of the high parton density in small x region a standard 
methods of perturbative QCD can not be used 0, ^ In this point the semihard 

approach [|[ |l^ differs from the works ^ |l^ where after choice of structure function 
at starting value Qq the QCD evolution is calculated using certain equations. 

In the region where the transverse mass of heavy quark M± = a/M^+p^ < < one 
need to take into account the dependence of the photon-gluon fusion cross section on the 
virtuality and polarizations of initial gluon. Thus in semihard approach the matrix elements 
of this subprocess differ from the ones of SPM (see for example Jt], [Tst). 

As for the SPM calculations of the next-to-leading (NLO) cross section of the photopro- 
duction of heavy quarks the review may be found in papers [|[ [T^ . The results of ||] have 
been confirmed by Smith and Van Neerven The further results for the electro- and 

photoproduction of heavy quarks are obtained in Refs. g 0. Authors of Ref. 
note that for beauty quark production the NLO corrections are large and various estimates 
of corrections lead to theoretical uncertainties of the order of a factor of 2 to 3. As it was 
estimated in ||2^ the total cross section for beauty quark production at HERA will be only 
a few tens of percent large than the one-loop results Q . 

For charm quark theoretical uncertainties are even higher. It is known that there is also 
strong mass dependence of the results of calculations for charm quark production. Since the 
mass of the charm quark is small for perturbative QCD calculations the resumme procedure 
is less reliable for charm quark production p2| . 

In Refs. we used the semihard approach to calculate the total and differen- 

tial cross sections of the heavy quarkonium, J/vP and T, photoproduction. We obtained 
the remarkable difference between the predictions of the semihard approach and the SPM 
especially for p±- and z-distributions of the J/^' mesons at HERA. 

In present paper we investigate the open heavy quark photoproduction processes in 
the semihard approach, which was used early in Ref. for calculation of heavy quark 
production rates at hadron colliders and for prediction of J/^ andT photoproduction cross 
sections at high energies in Refs. [^3[ p^ . 

Some of the results derived in this paper have already been presented in Ref. pq| . 
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2 QCD Cross Section for Heavy Quark Elec- 
troproduction 

We calculate the total and differential cross sections (the p±^ and rapidity distributions) 
of charm and beauty quark photoproduction via the photon-giuon fusion QCD subprocess 
(Fig.l) in the framework As for the SPM calculations of the next-to- leading (NLO) cross 
section of the semihard {k±^ factorization) approach 1^, 0. First of all we take into account 
the transverse momentum of gluon q2_L, its the virtuality ~ ^<72± the alignment of 
its polarization vectors along its transverse momentum such as = q2_L^/ \ 92i | 0, [T^ . 
Let us define Sudakov variables of the process ep — > QQX (Fig. 2): 

pi = aiPi + /3iP2 + Pi± P2 = ^2^1 + /32P2 + P2_L 

qi = XiPi + qii_ q2 = X2P2 + q2± (1) 

where 

22)1x2 22 22 

Pi and p2 are 4-momenta of the heavy quarks, qi is 4-momentum of the photon, q2 is 4- 
momentum of the gluon, pi±, P2±, qi±, 92_l are transverse 4-momenta of these ones. 
In the center of mass frame of colliding particles we can write Pi = {E,0,0, E), P2 = 
{E,0,0,-E), where E = ^/2, Pi = P^ = and (P1P2) = s/2. Sudakov variables are 
expressed as follows; 

Mi± , M2± . ^ 

ai = ^^exp(yi) a2 = exp(y2) 

Pi = ^^exp(-2/J /32 = ^^exp(-?/2): (2) 

where M{^2^ = + Pi 2_L^ Vi 2 rapidities of heavy quarks, M is heavy quark mass. 
From conservation laws we can easly obtain the following conditions: 

gii + 92^ Pi± +P2±, xi=ai+a2, X2=f3i+l32 (3) 
The differential cross section of heavy quark electroproduction has form: 

{ep QQX) = / dyldr, ^ ^ \M\'<^e{:ci C,!^)%{x2,qh) (4) 



(Ppi± J TT TT 167r^(a;ia;2s)^ 

For photoproduction process it reads: 



(Ppi± J TT 167r^(sa;2)^a2 

We use generalized gluon structure function of a proton ^p(x2,q2±) which is obtained in 
semihard approach. When integrated over transverse momentum q2± {q2± — (0, (7*2X70)) 
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of gluon up to some limit it becomes the usual structure function giving the gluon 
momentum fraction distribution at scale Q^: 



y <i>p(a;,g2\)rf9l± = a;Gp(a;,g2). (6) 



In our calculation we use the following phenomenological parametrizations | pj| : 

%{x,i^^) = - ^)Vp(^,92\), (7) 

gL) = ^1(1 - xf !^{x, g^), (8) 

where 

,2 ^ '2i 



\2 „2 ^ „2 

9| 



/p = (#^)^ 92\>'?^W, (9) 



and + exp(3.56v'ln(a;o/a;)), = 2GeV'^ , A = 56 MeV, xq = 1/3. The 

normalization factor C ~ 0.97 mb of the structure function ^p{x,q^^ in (7) was obtained 
in p^l l where &6-pair production at Tevatron energy was described. Thus the normalization 
constant C includes the so-called soft /C-factor, which takes into account as corrections due 
to soft gluon radiation. The typical value of ii'-factor in hard hadron-hadron interactions 
\% K — 2~ 2.5. As it was noted in Ref. |^ the choice of C = 0.97 mb corresponds to the 
upper limit of heavy quark hadroproduction cross sections. For the normalitation factor C\ 
in (8) we shall use value of C = 0.65 mb, which gives better description of existing data on 
charm quark photoproduction at fixed target energies \fs^^ = 10 — 30 GeV. 

In (4) <^e{xi,q\^^) is the well known virtual photon spectrum in Weizsacker- Williams 
approximation [E9 before the integration over g^j^: 



'^ixi,qi^)^—[ 10 

27r xiqfj^ qlj 



In this paper we cosider only heavy quark photoproduction processes. 

The effective gluon distributions xG{x,Q^) obtained from (7), (8) increase as x~'^" at 
not very small x (0.01 < x < 0.15), where uiq — 0.5 corresponds to the BFKL Pomeron 
smg ularity @. This rise goes continuously up to a; = xq, where xq being a solution of the 
equation qQ{xo) = Q^. In the region x < xq there is the saturation of the gluon distribution: 
2;G(2;,Q2) ~ CQ2 or CiQ2. 

The square of matrix element of partonic subprocess 7*5* — » QQ can be written as 
follows: 

|.,r|2 TR 2 2 r N2r 1 1 n , "2^1 , 0^1^23 n^\ 

\M\ ^16n e^aMx.-^s) [ (^ _ m^)(£- AP) ' + + ^ ^^^^ 
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For real photon and off-shell gluon it reads: 

|M|2 = 16^2g2 ( )2[_^ ^ ^(^L^ _ ^^^)% (12) 

where 02 = 1 — ai and s, u are usual Mandelstam variables of partonic subprocess 

s = (pi + P2f = {qi + q2f , i = ipi- qif = {P2- q2f, 

u = bi -92)' = (P2-'?1)', s + i + u = 2AP + ql^ + qi^. (13) 



3 Discussion of the Results 

The results of our calculations for the total cross sections of c- and 6-quark photoproduction 
are shown in Figs. 3 and 4. Solid curves correspond to the semihard approach predictions 
and dashed curves correspond to the SPM results with the GRV LO parametrization of the 
gluon distribution [Q. We used in our calculations rric — 1.5 GeV and = 4.75 GeV. 
Fig. 3a shows us that the solid curve for charm quark cross section obtained in semihard 
approach with parametrization (7) describes new ZEUS data as well as data from earlier 
fixed target experiments in the range 15 < %fa < 30 GeV. The curve of the parametrization 
(8) (Fig. 3b) describes better the experimental data at low energies (6 < ^/s < 30 GeV ) but 
gives worse description ZEUS data ||] . There are some differences between the theoretical 
curves and ZEUS experimental data (Figs. 3). If these ones will remain valid for new data, 
then its will show to essential resolved photon contribution ^ and a sensitivity of 
charm photoproduction cross section to the photon structure function [p9| . 

We would like to note that our results for the parametrization (7) coincide with the 
results from Ref. for pointlike component of the charm photoproduction cross section 
for MRSD'_ parton density at fi — rric and ones for parametrization (8) correspond to the 
charm photoproduction cross section in Ref. for MRS A parton density. 

At the energy range \/s^p > 500 GeV the saturation effects in the gluon distribution 
function are important for charm quark photoproduction in semihard approach (Fig. 4): 
the c-quark cross section grows more slowly at these energies (solid curve). The beauty 
quark photoproduction cross section predicted by semihard approach is larger than the one 
predicted by SPM at all energies (Fig. 4). Our conclusion for b-quark photoproduction 
is the same as in Ref. where b-quark hadroproduction was described using semihard 
approach. 

The p± distributions for c- and 6-quark photoproduction in the semihard approach (solid 
curves) and in the SPM (dashed curves) at the energy y/s^^ = 200 GeV are shown in Fig. 5. 
The curves are obtained in the semihard approach for charm quark photoproduction show 
the saturation effects in low p± region {pj_ < 2 Gev/c). At high p± region (3 < pj^ < 
20 GeV/c the heavy quark photoproduction p± distributions obtained in the semihard 
approach are higher than the ones of the SPM (with GRV LO parametrization of gluon 
distribution). This behaviour of c-quark p± distributions in the k± factorization approach 
results from the off mass shell subprocess cross section as well as the saturation effects of 
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the gluon structure function, because of ~ qo{x)^ at HERA energies. In the case of 
beauty quarks ml » qo{x)'^ and similar effects are absent. 

Fig. 6 show the comparison of heavy quark rapidity distribution (in the photon-proton 
center of mass frame) obtained in the different models: solid curve shows the y distributio 
semihard approach, dashed curve shows the one in SPM. The effects discussed above are 
sufficiently large near the kinematic boundaries, i.e. at big value of We see that 

the difference between solid and dashed curvers can't be degrade at all y* via change of 
normalization of both models. 

4 Conclusions 

We showed that the semihard approach describes experimental data for the charm pho- 
toproduction cross section at low and HERA energies, leads to the saturation effects for 
the total cross section of charm quark photoproduction at y/s^^ > 500 GeV as well as to 
exceeding over the SPM prediction for beauty quark total cross section at the energy range 
of \/s^p = 100 — 500 GeV, and predicts a marked difference for rapidity and transverse mo- 
mentum distributions of charm and beauty quark photoproduction, which can be studied 
already at HERA ep— collider. 
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Figure captions 

1. QCD diagrams for open heavy quark photoproduction subprocesses 

2. Diagram for heavy quark elecrtoproduction 

3. a) The total cross section for open charm quark photoproduction: solid curve - the 
semihard approach for the parametrization (7), dashed curve - the SPM. The solid 
dots are the ZEUS measurements |^ and the open dots are earlier measurements from 
fixed target experiments. 

b) The solid curve - the semihard approach for the parametrization (8). The experi- 
mental data as in Fig. 3a 

4. The total cross section for open charm and beauty quark photoproduction: curves as 
in Fig. 3a 

5. The p± distribution for charm and beauty quark photoproduction at ^s-^p = 200 
GeV: curves as in Fig. 3a 

6. The y* distribution for charm and beauty quark photoproduction: at y/s-yp = 200 
GeV: curves as in Fig. 3a 
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Fig. 1 




Fig. 2 
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